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5 ' Abstract 
< 

■ It is a challenging problem to construct an efficient quantum algorithm which can compute the 
Jones' polynomial for any knot or link obtained from platting or capping of a 2n-strand braid. 

, We recapitulate the construction of braid-group representations from vertex models. We present 

in 

the eigenbases and eigenvalues for the braiding generators and its usefulness in direct evaluation 
of Jones' polynomial. The calculation suggests that it is possible to associate a series of unitary 
operators for any braid word. Hence we propose a quantum algorithm using these unitary operators 
as quantum gates acting on a 2n qubit state. We show that the quantum computation gives Jones' 
polynomial for achiral knots and links. 
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I. INTRODUCTION 



Quantum algorithms have proven to be more efficient than classical algorithms in solving 
a number of problems Q|. For instance, Shor's quantum algorithm can factor numbers 
exponentially faster than classical algorithms. Grover's search algorithm does not change 
the complexity class, but provides significant speed up for large databases. We would like 
to explore the power of quantum algorithms in the context of knot theory. 

Classification of knots and links in a three-dimensional space is one of the open problems. 
Jones introduced a recursive procedure for determining a polynomial relation for these knots 
and links. Jones' polynomials do classify some knots and links Q]. There are other general- 
ized polynomials which improve the classification but none of them have achieved complete 
classification |3L It is known that the evaluation of Jones' polynomial classically is a j^V 
hard problem Hence it will be interesting to study the computation of knot and link 
polynomials using a quantum algorithm. 

There are diverse approaches in physics to obtain polynomials for knots and links. Fol- 
lowing Alexander's theorem, any knot can be viewed as a closure or capping of an n-strand 
braid. Therefore the polynomials for knots and links can be determined by studying repre- 
sentation theory of braid groups B n . The common ingredient in these approaches is to find 
different representations of braid groups B n . We now present a brief summary of some of 
these approaches: 

1) N-state vertex models which are two-dimensional statistical mechanical models 
where the bonds of the square lattice carry spin n/2 representations of 377(2). The number 
of possible states of spin n/2 is denoted by N — n + 1. The properties of these models are 
described by the so-called i?-matrix which is an N 2 x N 2 matrix. The number of nonzero 
elements in the i?-matrix for iV-state vertex models is given by m = ^S 1 1 2s 2 + N 2 . In 
the literature, the vertex models are referred to either as iV-state vertex model or m-vertex 
models; both are equivalent. For example, 2-state vertex models carry spin 1/2 on the 
lattice bonds and they are equivalently called as six-vertex models where "six" denotes the 
number of non-zero R- matrix elements. In Ref. braid group representations and knot 
polynomials from the i?-matrices of N = 2, 3, 4-state vertex models were obtained. 

2) Chern-Simons gauge theory is a topological field theory which provides a natural 
framework for the study of knots and links p. The knot polynomials are given by the expec- 
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tation value of Wilson loop observables. In particular, the Jones' polynomial corresponds to 
the Wilson loop carrying spin 1/2 representation in SU(2) Chern-Simons theory. Clearly, ar- 
)itrary representations of any compact gauge group G can result in generalized polynomials 
3]. The polynomials are in the variable q which is a function of the coupling constant k and 
the rank of the gauge group. The field theoretic polynomials were obtained by exploiting 
the connection between Chern-Simons theory on a three-manifold with boundary and the 
corresponding Wess-Zumino-Witten conformal field theory (WZW) on the boundary. The 
polynomials crucially depended on various representations of the monodromy or braiding 
matrices in the WZW models. 

Recently, Freedman et al. [9] have attempted simulation of topological field theories by 
quantum computers. The topological quantum computation proposed in Refs. j^-jlll is at 
a mathematically abstract level. It exploits the connection between fractional quantum Hall 



bj. In Ref. |rj, the 



states and Chern-Simons theory at the appropriate integer coup 

3) State sum method of obtaining bracket polynomials 
construction of a unitary representation has been shown for the three-strand braid B3. 
Further within this approach, it has been shown that it is not possible for a quantum 
computer to evaluate the knot polynomial. However, for a specific choice of the polynomial 
variable, the linking number can be determined 

Our aim is to determine the Jones' polynomial for any knot or link obtained from braids 
using a quantum algorithm. For this purpose, we need to determine matrix representation 
for braid generators. We will recapitulate the construction of braid group representation 
from the six-vertex model. Then, we can determine the eigenvalues and eigenstates for 
these braid generators. This exercise suggests that we can associate product of unitary 
operators for any braid word. Hence the unitary transformations, corresponding to any 
braid word, can be implemented using quantum gates. It is important to stress that the 
quantum computation in this paper is crucially dependent on the mapping of any braid word 
to a product of unitary operators. 

The polynomials for knots and links can be directly computed by choosing a suitable 
eigenbasis of the braiding matrices. Essentially, a quantum algorithm will determine the 
probability of finding unitarily evolved initial state in a final state (\out)): 

P K = \( out \U B \in )| 2 (1) 
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where Ub represents the series of unitary operators corresponding to the braid word. For 
the braiding matrices obtained from the six-vertex model, the above matrix element gives 
the modulus-square of the Jones' polynomial (up to an overall constant). 

This paper is organized as follows. In Section |H1 we present a general method to find 
representation of braid groups using N-state vertex models. We discuss in detail the six- 
vertex model, braiding eigenvalues and eigenstates. Using these eigenstates, we evaluate 
Jones' polynomial. In Section IIII1 we present a method to perform the evaluation of the 
modulus-square of the Jones' polynomial as a quantum computation by considering any 
knot or link as a composition of cups, a series of braiding operations and caps. In the 
concluding section, we summarize the results obtained and discuss the significance of the 
quantum algorithm. 

II. N-STATE VERTEX MODEL 

In this section, we review the construction of braid group representations from N-state 
vertex models In order to compare the eigenstates of the braiding operator with the 
qubit states, the six- vertex model (spin 1/2 on the bonds of the square lattice) is relevant. 
Hence, we shall present the explicit form of the i?-matrix and the braid matrix for the 
six-vertex model. 

As mentioned in the introduction, vertex models are two-dimensional statistical me- 
chanical models with the spins n/2 lying on the bonds of a square lattice. The proper- 
ties of these models are described by the R- matrix elements between edge states (mi,m 2 ) 
and (711,712): Rm™m 2 ( u ) where u is the spectral parameter. Here 7711,7712,711,712 take values 
n/2, n/2 — 1, ... — n/2. The integrability condition of these models requires the following 
equations to be satisfied: 



where Xi(u)s are called Yang-Baxter operators and the relation (J2J) is called Yang-Baxter 
equation. The explicit form of Xi(u) in terms of the .R-matrix elements is given by 





(2) 
(3) 




(4) 
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where 1^) i s the identity acting at the j-th position and E is a matrix such that (E nim2 ) pq 
8n lP 8m 2 q- The solution to Eq. (J2J) can be written in a compact form 



i?™ 1 ™ 2 Cm) 



J 



j=o ™2 m 



Aj(tt) ' 1 



J 



777,2 mi m 



(5) 



where the terms in parenthesis are the quantum Clebsch-Gordan coefficients (q-CG) 
which are nonzero if and only if m takes a value in the range m — — J, — J + 1, . . . + J and 
satisfies the condition mi + 777.2 = m = n\ + 77,2- Here Xj(u) is given by 

jl+j'2 J 



I j sinh(s/i + u) ] [ sinh(/yU — w) 



(6) 



S=J+1 



where the q-CG coefficient variable q = e 2fl . These solutions are spectral parameter depen- 
dent solutions. The explicit form of the i?-matrix for the six-vertex model is 



n x n 2 
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We are interested in the construction of braid group representations from these Yang- 
Baxter operators Xi(u) |5[. The defining relations of the braid group generators are exactly 
same as Eqs. (|2I3|) in the absence of the spectral parameter. In other words, the braid group 
generators hi e B n (spectral parameter independent operators) are obtained from Xi(u) by 
taking the spectral parameter u — > 00: 



lim Xiiu) 

u^oo 



(7) 



Let us define a quantity a™ 1 ™^ as follows: 



lim P nin2 = (T n i n 2 



(8) 



The explicit form of the a^^-matrix elements in this limit turns out to be: 
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The above matrix elements indicate that we can choose a basis for the n-strand braid 
as |mxm2 . . . mjm i+1 . . . m n ) where m^s take values ±| which are the two possible quantum 
states of the spins s^'s = |. Hence rrii denotes a single qubit basis and the above basis state 
is an n-qubit basis state. From Eq. (J3J, we can write the action of the braiding operator hi 
on such an n-qubit state: 

b i \m 1 m 2 . . . mim i+l . . . m n ) = ^ Oji l™i™2 • • • • • • m n ) , (9) 

where n^n^+i take values ±|. Clearly, this equation is not an eigenvalue equation of the 
braiding operator. We need to diagonalise the cr^^n 2 matrix to determine the eigenbasis of 
the braiding operator 6j. As we did in the case of Eq. (jSJ), we can diagonalize the spectral 
parameter independent cr^ 1 ^ using the quantum Clebsch-Gordan coefficients matrix. The 
corresponding eigenstates will be similar to the coupled states obtained from uncoupled 
states using quantum CG-coefficient matrix jis| : 

| Jim) = ^ | * +1 J ||mi,m i+ i) , (10) 

mj,m i+ i m i m i+l vr ^ 

where Si, s$+i are the spin | states and Jj 6 Sj ® s i+ i = 0, 1. Thus we obtain a 2-qubit state 
\Jiw) which may be an entangled state. 

The form of the diagonal a matrix in the coupled basis is: 



\Jm) 


111) 


|10) 


11 -1) 


|0 0) 


\Jm) 
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where Ai = 1 and Ao = —q. These eigenvalues are equal up to an overall normalization to 
the eigenvalues of the Wess-Zumino-Witten model monodromy matrices. 

We observe that the eigenvalues of the a matrix on the coupled states | Jm) depend only 
on J and not on m. Therefore, we can suppress the m dependence on the eigenstates of the 
braiding operator and equivalently write it as a tensor product state involving the spin | 
placed on the bonds of the six-vertex model. That is, 

\JiTn) = \Ji) = \(si <g> s i+1 )j t } . (11) 

Even though we have explicitly diagonalized the a matrix, we must remember that all the 
braid group generators b^s cannot be simultaneously diagonalized. The spectral parameter 
independent form of Eqs. (jZEl are the defining relations of the braid group B n which implies 
that we can simultaneously diagonalize either 6 2 j's or 6 2 i+i 's. 



A. Representation Theory of Braid Groups in Knot Theory 

In this subsection, we would like to address the eigenvectors and eigenvalues of braid 
generators from the viewpoint of obtaining polynomial invariants of knots from platting or 
capping of braids. 

It is well known that knots from braids are not unique. That is, braids related by Markov 
moves I and II give rise to the same knot. These two moves indeed completely remove the 
non-uniqueness. So the construction of polynomial invariants for knots must be such that the 
polynomial does not change under Markov moves. One such procedure for knots obtained 
from closure of braids has been presented in . 

We will use the eigenstates of the braiding operators to directly compute the polynomial 
invariant for any knot obtained from closure or capping of an n-strand braid. In order to 
remove the non-uniqueness due to Markov moves, we place orientations on the strands of 
the braid. Further, we introduce a correction factor to the braid eigenvalues obtained from 
six-vertex model such that the polynomial does not change under Markov moves I and II. 
The correction factor on braiding eigenvalues depends on the relative orientations between 
the two strands. For right-handed half-twists between strands of parallel orientation, the 
braiding eigenvalues Aj + ' ) are 

A™ = -q^ ; A« = gV2 . (12) 
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Similarly, for right-handed half twists between strands of antiparallel orientation, the braid- 
ing eigenvalues Xj ■* are 

A<f } = 1 ; AjfJ = -q- 1 . (13) 

The eigenvalues for left-handed half-twists are inverse of the right-handed half-twists eigen- 
values. 

Suppose we consider any knot obtained from capping of a 2n-strand oriented braid. 
Clearly, capping is possible if the number of outgoing strands is equal the number of incoming 
strands in the oriented braid. In other words, the quantum states m^'s on the strands (JHJ) 
should be such that 

2n 

Y,rrii = 0. (14) 
i=i 

Therefore to study knots from braids, only the subspace of the states |mim 2 ■ ■ - Tn^n) satis- 
fying Eq. (|14jl needs to be considered. Hence the construction of the eigenstates of braiding 
operators should be consistent with Eq. p4jl . 

We will now present the eigenstates of the braiding matrices which will enable direct 
evaluation of knot polynomials. For a 2n-strand oriented braid, we can write the most 
general eigenbasis of braiding operators 6 2 i+i with eigenvalue Xj 2i+1 , for all z's, as: 

\<f>({J2i+lUk})) H( •••(•••( (si ® s 2 )ji ® (s 3 ® «4)j 3 )h ^ 
<8> . . . (S 2 i+1 ® S 2 i+ 2 )j 2!+1 )i I ® • • • («2n-l ® S 2 „)j 2n _ 1 ) ) 

Recall that the appropriate braiding eigenvalues (J12I13)) need to be substituted depending 
on the relative orientations of the two strands involved in braiding and the handedness. The 
brackets within the basis kets should be identified with the notation in Eq. 1)11)1 ■ That is, 
|(si ® s 2 )j 1 ) = I Ji), |((si <g) s 2 )j 1 ® (s 3 g) s 4 )j 3 ) il ) = I (Ji (g) ^3)^1), and so on. Note that the 
final combined state in the above basis is chosen to be spin which is essential to satisfy 
the condition (JH)) to describe knots from closure or capping of braids. 

In the similar fashion, we can write a different eigenbasis for braiding operators 6 2 j with 
eigenvalue Xj 2i for all i's: 

\k{J2iUn})) =!((••• ( • • • ( ( si <g> (s 2 ® s 3 ))j 2 ) ri <g> . . . 

(16) 

(S 2 i (g> S2i+l)j 2i )n-l ® • • • ( s 2n-2 <g> S 2n _l ) J 2n _ 2 ) r „_ 2 (g) S 2n ) ) 

In order to achieve the final spin state, we require r n _ 2 = |. Incidentally, these two bases 
are equivalent to the conformal blocks in Wess-Witten conformal field theory. 
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The two different bases ()15ll6j) are related by an orthogonal (unitary) duality matrix 



JjK+i},{Ji})) = a ({^ l +i},{i l }),({J2 l },{n}) IfewMn})) ■ 



(17) 



The duality matrix can be written in terms of products of SU(2) q quantum- Racah coefficient 
matrices |l7|: 



a({J 2i+ i},{h}),aJ2 t },{n}) 



n-2 

tl,...t n -2 1=1 



n-2 

n a *f-a 

m=0 



ti-1 ^2i 
r i $2m 



h-1 s 2i+l 
S2i+2 U 

t m S2m+2 
s 2m+3 t m+ i 

where the closed form expression for the quantum-Racah coefficient matrix is Q| 



n-2 
m=0 



(18) 



si s 2 

S3 S 4 



\S1+S2+S3+Si 



y/[2j + l][2l + l]A( ai , s 2 ,j)A(s 3 , s A ,j)A( 8l , s 4 , l)A(s 2 , s 3 , 1) 



x l) m [m + l]!{[m — si — s 2 — j]l[m — s 3 — s 4 — j]\[m — s± — s 4 — l]\ 

m>0 

x[m - s 2 - s 3 - l]\[si + s 2 + s 3 + s 4 - m]![si + s 3 + j + / - m]\ 

x[s 2 + s 4 +j + /-m]!}" 1 (19) 



where A(a, 6, c) = 
g-number defined as 



-a+6+c] ! [a-6+c] ! [a+b-c] ! 
[a+fe+c+1]! 



and the number in square brackets represents the 



q i/2 Z q -i/2 ■ 



With this detailed background on the representation theory of braid groups, we are now 
in a position to evaluate Jones' polynomial for any knot or link. 

B. Evaluation of Jones' polynomial 

Consider any knot or link as shown in Fig. 1 which is technically called platting or capping 
of braids. The knot type is determined by the braiding region denoted as the shaded region 
in Fig. 1. It involves a sequence of braiding operations which is usually written as a braid 
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B = Braid Word 
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FIG. 1: Viewing Knots as a Capping of Braiding Sequences 

word with appropriate orientations. Taking the initial state and final state as |0({o},{o}))j 
the following matrix element 

Vk = (0({o},{o})|B|0( {O },{o})) (20) 

gives the Jones' polynomial (up to an overall normalisation). To see the polynomial form 
in variable q, we can substitute the eigenvalues of the braiding generators and the unitary 
duality matrices elaborated in the previous subsection. 

For the variable q being root of unity, we can equivalently denote the sequence of braiding 
operations B e £> 2n as a product of unitary 2 2n x 2 2n matrices. We will now illustrate this 
mapping for two examples involving four-strand and six-strand braids. As shown in Fig. 2 (a), 
the braid word 

B a = blbfbfbl e B, , 

where bi denotes braiding between anti-parallely oriented i-th. strand and i+l-th strand and 
bj denotes braiding between parallely oriented j-th strand and j+l-th strand. Following the 
rules of the representation theory of braids, the braid- word in Fig. 2(a) will be mapped to 
a product of 2 4 x 2 4 unitary matrices as follows: 

B a = XJ Ba = ({a f (b 2 ) a f g(bi, b 3 ) a h(b 2 ) at} 2x2 © Ii4xi4) 24x24 (21) 

where a is the duality matrix and the matrix elements of the matrices f (b 2 ), g(b 1; b 3 ) and 
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FIG. 2: Two examples (a) B a € B 4 , (b) B b e B& 

h(b 2 ), for the braid word B a , will be 
f(ba) JlJa = WA^) 3 , g(b 1 ,b 3 ) Jl;j2 = S JlJ2 (\^)- 4 , h(b 2 ) JiJa = 6 JlJ2 (\^r . (22) 

The entries in these diagonal matrices are functions of the braiding eigenvalues. The relative 
orientations and the number of crossings between the respective braids are incorporated in 
writing the functional form. If q is a root of unity, these diagonal matrices are unitary. Thus 
we see that the braid word B a e £> 4 can be equivalently represented as a product of unitary 
2 4 x 2 4 matrices. Similarly, the braid word B\, corresponding to Fig. 2(b) will be 

B b = (6j 1 6 4 )(Sr 2 63 3 65 2 )(S 2 S2)6 1 & 2 e B 6 . 

The equivalent description involving product of 2 6 x 2 6 unitary matrices will be 

XJ Bb = [(a f (b 2 , b 4 ) a f g(b x , b 3 , b 5 ) a h(b 2 , b 4 ) a f fi(bi) a gi(b 2 ) a f ) 5x5 © l 5 9x59]26 X 26 

(23) 
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where the duality matrix elements are denoted as a(j 1 ,j3,j 5 =z 1 ),(j2,j4,ri) and the remaining 
5x5 diagonal matrices for the braid word Bb will be 

f(b 2 ,b 4 ) (J , ,J/,ri),(jr a ,J 4 , ri ) = ^J 2 J2^J 4 ,J4^,r 1 (Aj2 ) ) _1 A^ ) 
g(bi,b 3 ,b 5 ) (J , iJ , iJ , ))(Ji)j3jj6) = ^j ^^.jSX, ) 2 (A«) "(X,;) 2 
h(b 2 ,b 4 ) (J , iJ , r , )i(j2iJ4iri) = 5^ j2 5j kM 5 r>1 \ { ;?{\ { ^) 2 

f l( b l)(Ji,J5,^),(Jl,Js,J6) = 5 J^i 5 4^3 5 4,J s Ajt 

gl( b a)(J i) Ji,ri),(J a ,J4,ri) = ^ 2 J2^J 4 ,J4^,r 1 A^ ) . (24) 

In a similar fashion, we can find the unitary representation Ub for any braid word 
B G Bm in terms of products of 2 2n x 2 2n unitary matrices involving duality matrices 
[ a ({j 2 i + i},{ii}).({-fe},{M) ® l]2 2n x2 2n an d diagonal braiding matrices. These unitary represen- 
tations play the role of quantum gates in the quantum computation of Jones' polynomial. 



III. QUANTUM COMPUTATION 

In this section, we attempt to compute the Jones' polynomials, for knots and links ob- 
tained from platting or capping of In strand braid as shown in Fig. 1, through a quantum 
algorithm. We have already elaborated in the previous section that we can associate Ub 
(product of unitary matrices) for every braid word B G £>2 n - The quantum algorithm in- 
volves the following steps: 

Step 1: Let the initial 2n-qubit state be |0) (\cup)). 

Step 2: We perform the sequence of unitary operations corresponding to the braid word 
B G i?2n in Fig. 1. The unitarily transformed state will be 

|*) = U B |0). (25) 

Step 3: Finally, we determine the probablity of the unitarily evolved state |\&) in a specific 
final state |/) as 

I(/|Ub|0)| 2 . (26) 

Taking the final state to be \cap) = |0), we get the modulus square of the Jones' polynomial 
(up to an overall normalisation)^ (J2*UJ| . 
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For a subclass of knots(links) called achiral knots(links), Vk is unchanged under q — > q^ 1 . 
In other words, the matrix element (0|£/b|0) will be real. For these achiral knots and links, 
the quantum algorithm directly gives the Jones' polynomial (up to an overall normalisation). 

IV. SUMMARY AND DISCUSSIONS 

In this paper, we have presented matrix representations for braiding matrices from six- 
vertex models. We have discussed the representation theory of braids, namely, the eigenbasis 
and eigenvalues of the braid generators obtained from six- vertex models. The explicit evalu- 
ation of Jones' polynomial, for any knot/link from braids, is presented. From the evaluation, 
we have shown that we can associate a series of unitary operators for any braid word. This 
is the significant result of the paper enabling quantum computation. We have demonstrated 
a quantum algorithm, involving these unitary operators, which can determine the modulus 
square of the Jones' polynomial for any knot or link. The algorithm gives Jones' polynomial 
for achiral knots and links. 

We must realize that the quantum computation essentially determines the probablity 
of unitarily evolved initial state in a specific final state. Further, the number of unitary 
operators is dependent on the braid word and at most equal to twice the length of the braid 
word. 
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